for i in [2, 2,4,8,16,32,64,128,256,512,1024,2048,4096]
    T1 = time()
    global nels = i
    include("main1.jl")
    println("neq = ", neq, " cputime = ", time() - T1)
end
